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ABSTRACT 


This  thesis  presents  the  bandpass  inverse  fast-Fourier  transform  (IFFT)  filter  bank 
and  the  multirate  digital  filter  bank  techniques  to  synthesize  test  point  waveforms  from 
constituent  waveforms  recorded  by  two  instruments  as  part  of  an  aircraft  electromagnetic 
hardness  evaluation  test.  The  component  waveforms  are  recorded  by  two  separate 
measurement  systems  (High-Powered  Pulse  Waveform  (HPW)  in  the  time  domain  and 
Continuous  Sweep  Waveform  (CSW)  in  the  frequency  domain)  imder  two  different  aircraft 
orientations  (parallel  and  perpendicular).  Data  from  two  orientations  are  combined  using 
the  sinusoidal  modeling  algorithm  (SMA).  The  tree-structured  filter  bank  with  power 
symmetric  overlap  method  and  the  bandpass  IFFT  with  spectral  concatenation  method  are 
developed  to  further  combine  these  waveforms  with  an  overlapping  frequency  spectrum  to 
produce  the  corresponding  synthesized  test  point  waveform. 
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I.  INTRODUCTION 


A.  WIDEBAND  TRANSIENT  WAVEFORMS 

Testing  of  aircraft  for  electromagnetic  hardness  evaluation  involves  measurement  of 
wideband  electromagnetic  transient  waveforms.  It  is  difficult  to  devise  a  single 
measurement  system  for  such  wideband  signals.  As  a  result,  wideband  measurements  are 
often  performed  by  more  than  one  instrument.  These  separate  measurements  covering 
different  frequency  bands  are  then  combined  together  to  synthesize  an  equivalent  single 
waveform  in  both  time  and  frequency  domains.  The  data  sets  used  in  this  thesis  are 
representative  of  measmements  obtained  from  an  aircraft  hardness  test.  Measurements  are 
typically  made  at  specific  test  points  inside  the  aircraft.  Four  waveforms  are  collected  at 
each  test  point:  two  waveforms  from  the  High-Powered  Pulse  Waveform  (HPW) 
measurements  and  two  waveforms  from  the  Continuous  Sweep  Waveform  (CSW) 
measurements.  Each  of  these  waveforms  is  obtained  from  either  parallel  or  perpendicular 
orientation  of  the  aircraft.  The  HPW  data  set  is  recorded  in  the  time  domain  while  the 
CSW  is  measured  in  the  frequency  domain.  The  two  sets  of  measurements  cover  different 
frequency  bands  with  a  small  overly  in  the  80-100  MHz  range.  Different  sampling  rates 
are  used  to  record  the  data,  and  the  frequency  domain  measurements  are  actually  recorded 
in  six  different  frequency  bands  at  varying  bin  widths. 

B.  GOAL  OF  THESIS 

The  goal  of  this  thesis  is  to  develop  techniques  for  combining  component 
waveforms  measured  at  a  test  point  in  the  aircraft  to  produce  a  synthesized  test  point 
waveform.  These  waveforms  are  measured  imder  two  different  orientations  and  over  two 
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different  frequency  bands.  In  previous  work,  some  algorithms  were  designed  by 
Winnenberg  [Ref.  1]  and  Martin  [Ref.  2]  to  produce  synthetic  signals  with  close  statistical 
similarity  to  original  signal.  In  this  thesis,  the  bandpass  IFFT  and  the  multirate  digital  filter 
banks  are  employed  to  combine  the  measured  data  to  produce  a  synthesized  test  point 
transient  waveform  which  still  retains  the  basic  properties  of  the  component  signals.  There 
are  four  waveforms  to  be  combined:  two  HPW  and  two  CSW  waveforms,  which  together 
encompass  the  spectrum  from  0.5  MHz  to  1  GHz. 

The  concept  of  wideband  analysis  as  it  applies  to  digital  signal  processing  (DSP) 
can  be  explained  as  a  method  of  dividing  a  large  frequency  range  into  smaller,  more 
manageable  subbands  [Ref.  3].  Thus,  the  analysis  will  essentially  divide  the  original 
signal  into  groups  of  signals  which  retain  the  frequency  characteristics  of  the  original 
signal.  FFT  and  multirate  filter  banks  are  used  to  perform  wideband  signal  analysis  and 
synthesis.  Since  a  method  of  physically  measuring  the  test  point  signal  covering  the  entire 
spectrum  of  interest  does  not  exist,  it  is  not  possible  to  directly  determine  the  accuracy  or 
effectiveness  of  the  work  presented.  Therefore,  the  norm  attributes  of  each  signal  are 
computed  to  determine  the  closeness  of  the  synthesized  test  point  signal  to  its  component 
waveforms. 

C.  THESIS  OUTLINE 

The  remainder  of  this  thesis  is  organized  as  follows.  Chapter  n  introduces  the 
properties  of  bandpass  signals  and  the  bandpass  IFF!  technique.  Chapter  IQ  discusses  the 
theory  and  design  of  multirate  filter  banks.  Chapter  IV  evaluates  the  performance  of  the 
bandpass  IFFT  and  the  multirate  filter  bank  for  a  given  set  of  test  point  measurements.  The 
synthesized  test  point  waveforms  are  evaluated  based  on  visual  observations  and  numerical 
comparison  iwing  norm  attributes.  Chapter  V  presents  conclusions. 
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II.  BANDPASS  SIGNALS 


A.  THEORETICAL  BACKGROUND 


It  is  essential  to  understand  how  to  efficiently  change  the  sampling  rate  of  bandpass 
(rather  than  lowpass)  signals  because  this  work  involves  M -channel  filter  banks,  which 
provide  the  ability  to  decompose  the  signal  into  M  subband  components. 

1.  Complex  Lowpass  Representation  of  Bandpass  Signals 

The  spectrum  of  bandpass  signals  is  nonzero  only  within  a  region  of  around 
some  carrier  (center)  frequency,  f^>W,  where  2W  is  the  bandpass  signal  bandwidth. 
Typically  a  bandpass  signal  is  narrow  band  if  the  bandwidth  of  the  signal  is  smaller  than  its 
midband  frequency,  i.e.,  f^>2W.  A  bandpass  signal  is  represented  in  general  as 

x{t)  =  A(r)cos[23t/cr  +  <l)(0]  (2.1) 

where  A{t)  is  the  envelope  of  the  signal,  and  (j>(?)  is  the  phase  of  the  signal  [Ref.  4], 

The  amplitude  spectrum,  |X(/)| ,  of  a  typical  bandpass  signal,  x{t) ,  is  shown  in 
Figure  2.1  (a)w  The  pre-envelope  of  the  bandpass  signal  is  defined  by 

x+(r)  =  xit)  +  jx{t)  (2.2) 

where  x{t)  is  the  Hilbert  transform  of  x{t) . 
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Figure  2.1  Complex  Lowpass  Representation  of  Bandpass  Signals.  After  [Ref.  6]. 


The  pre-envelope,  is  a  complex- valued  function  of  time  with  the  original 
signal,  x(t) ,  as  the  real  part  and  its  Hilbert  transform,  5c(t) ,  as  the  imaginary  part.  If 
X^if)  is  defined  as  the  Fourier  transform  of  and  X(f)  is  defined  as  the  Fourier 
transform  of  x{t) ,  then 


X+(/)  =  X(/) +  ;!(/) 


(2.3) 


From  Haykin  [Ref.  6],  it  is  known  that  X{f)  =  -jfsgn(/)X(/),  where  sgn(/)  is  the 
signum  function.  Then,  Equation  2.3  may  be  written  as 

Xt(/)  =  X(/)  +  sgn(/)X(/)  (2.4) 


Using  the  definition  of  the  signum  function.  Equation  2.4  becomes 

xHf)  = . 


mf), 

X(0), 

0. 


/>o 

/  =  0  (2.5) 

/<0 

where  X(0)  is  the  zero  frequency  value  of  X(f) .  From  Equation  2.5,  it  is  clear  that  there 
is  no  negative  frequency  component  in  the  pre-envelope  of  a  Fourier  transformable  signal 
as  shown  Figure  2.1(b). 

From  the  frequency  shifting  property  of  the  Fourier  transform,  the  pre-envelope 
may  be  expressed  as 


;ct(r)  =  (2.6) 

where  ,x(r)  is  a  complex- valued  lowpass  signal  whose  amplitude  spectrum  is  shown  in 
Figure  2.1(c). 

The  complex  envelope,  x{t) ,  of  a  given  bandpass  signal,  x{t) ,  may  be  determined 

as  illustrated  in  Figure  2.  l(a)-(c)  in  the  following  manner: 

Step  1.  Retain  die  positive  frequency  half  of  X{f)  centered  at  . 

Step  2.  Shift  it  to  &e  left  by  . 

Step  3.  Scale  it  by  a  factor  of  two. 

In  this  approach,  the  Hilbert  transform,  x{t) ,  does  not  need  to  be  computed.  The 
real  part  of  A:^’(r)  is  J[:(r)  =  Re{x{t)J  ] .  From  Equation  2.6,  any  real  bandpass 
signal,  x(t),  can  be  represented  in  terms  of  an  equivalent  complex- valued  lowpass  signal, 
x(t). 
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2.  Sampling  of  Bandpass  Signal 


It  is  known  that  if  a  lowpass  signal  is  sampled  at  twice  its  highest  frequency,  the 
original  signal  can  be  perfectly  reconstructed.  However,  it  is  possible  to  sample  a 
bandpass  signal  at  a  minimum  rate  of  twice  its  bandwidth  without  any  loss  of  information. 
Figure  2.2  describes  sampling  of  a  bandpass  signal. 


Figure  2.2  Frequency  Translation  of  Bandpass  Signals;  (a)  Spectrum  of  Bandpass  Signal 
and  (b)  Its  Lowpass  Representation.  After  [Ref.  7]. 

Consider  a  banc^ass  signal,  Xgpif) ,  which  contains  spectral  components  only  in 
the  frequency  range,  //  <  |/|  <  //  +  /^  •  denote  the  component  of  Xpp(f) 

associated  with  /  >  0  and  Xp(/)  denote  the  component  of  Xpp{f)  associated  with  /  <  0 , 
as  shown  in  Figure  2.2(a). 

Figure  2.2(b)  illustrates  that  by  lowpass  translating  (modulating)  X„(/)  to  the 
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band  0  to  and  Xp(/)  to  the  band  - to  0,  the  new  signal,  X^pif) ,  can  be  generated 
which  is  “equivalent”  to  Xgp(f)  in  the  sense  that  Xppif)  can  be  uniquely  reconstructed 
from  Xj^pif)  by  the  inverse  process  of  bandpass  translation.  However,  by  applying  the 
concept  of  lowpass  filtering,  it  can  be  shown  that  the  sampling  firequency  required  to 
represent  this  is  now  f^ip  ^  2/^ ,  which  can  be  much  less  than  f^pp  k  2(//  +  /^) . 

In  practice,  there  are  many  ways  in  which  the  combination  of  frequency  translation 
and  sampling  described  above  can  be  performed.  Several  modulation  techniques 
including  integer-band  sampling,  quadrature  modulation,  and  single  sideband  modulation 
are  desoibed  in  [Ref.  7]. 

B.  INVERSE-FOURIER  TRANSFORM  OF  BANDPASS  SIGNALS 


Each  subband  signal  of  the  CSW  waveform  is  recorded  at  a  different  sampling  rate 
with  no  measurements  for  negative  frequencies.  This  signal  may  be  considered  a  pre¬ 
envelope  of  a  bandpass  signal  which  is  Fourier  transformable.  Let  s(t)  be  one  of  the 
subband  signals  of  the  CSW  waveform,  which  is  recorded  in  the  frequency  domain. 

Because  there  is  no  data  available  on  the  negative  frequency  side,  we  handle  this  signal  as 
s'+Ct) ,  which  is  a  pre-envelope  of  5(t) .  The  amplitude  spectrum,  l5f(/)| ,  of  this  signal  is 
depicted  in  Figure  2.3(a). 

From  the  inverse  processes  of  Equations  2.2,  2.3,  and  2.6,  we  may  obtain  the 
subband  signal  in  the  frequency  domain  as  depicted  in  Figures  2.3(a),  (b)  and  (c)  in  the 
following  manner: 

Step  1.  Scale  the  pre-envelope  signal,  |5+  (/)| ,  by  a  factor  of  one  half. 

Step  2.  Obtain  its  mirror  image  and  locate  it  on  Ihe  negative  frequency  side 
centered  at  . 

Step  3.  Zero-pad  the  points  between  the  original  and  the  mirror  image  spectrum. 
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Figure  2.3  Transformation  of  a  Bandpass  Signal  from  the  Frequency-Domain  to  the  Time' 

Domain. 


In  Step  1,  care  must  be  taken  to  ensure  that  the  frequency  domain  data  satisfy  the 
symmetry  property  about  the  origin,  so  it  is  possible  to  recover  the  real  valued  time  domain 
data.  Figure  2.3(c)  shows  the  spectra  after  zero-padding.  This  signal  can  be  inverse 
Fourier  transformed  to  obtain  the  corresponding  time  domain  signal  for  further  analysis. 
These  steps  can  be  interpreted  as  a  bandpass  IFFT  operation  (see  Figure  2.3(d)). 

C.  BIN  WIDTH  CONVERSION  OF  BANDPASS  SIGNALS 


Figure  2.4  shows  a  typical  CSW  waveform,  which  is  recorded  in  the  frequency 
domain.  As  was  mentioned  earlier,  different  sampling  rates  were  used  in  each  subband  as 
shown  in  Figure  2.5.  Each  subband  signal  is  needed  to  be  interpolated  and/or  decimated 
to  realize  the  overall  CSW  waveform  with  a  constant  sampling  rate. 


Figure  2.4  Magnitude  and  Phase  Measurements  of  a  CSW  Waveform  (LAIOlOOl). 
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Let  A/j ,  A/2 . .  A/g ,  be  the  frequency  bin  widths  of  subband  signals.  Si,  S2,  ■ 
. . ,  Sg ,  respectively,  as  shown  in  Figure  2.5.  Table  2.1  shows  the  frequency  ranges  and  bin 
widths  of  the  six  subband  signals  of  the  CSW  waveform. 


5i  S2  S2  S4,  5g 


A/i 

A/2 

A/s 

A/4 

A/s 

A/6 

0  0.51  2.192  40.55  121.6  156.9  243.1  1000  MHz 

Figure  2.5  CSW  Subbands 


Table  2.1:  Subband  Frequency  Ranges  and  Bin  Widths  of  Test  Data 


Subband 

Signal 

Frequency  Band 
(Hz) 

Bin  Width 
(Hz) 

Si 

0.510e6  -  2.192e6 

0.0847e5 

S2 

2.192e6  -  40.55e6 

0.365  le5 

S3 

40.55e6-  121.6e6 

6.7375e5 

S4 

121.6e6-159.6e6 

19.500e5 

S5 

159.6e6  -  243.  Ie6 

26.060e5 

S6 

243.  Ie6  -  1000e6 

40.476e5 

Let  Ni,N2,..  .  ,Nq  be  the  number  of  sampled  data  points  of  subband  signals  5^ , 
S2, .  • .  .iSg,  respectively.  Subband  signals  52 -Sg  are  interpolated  to  achieve  a  uniform 
bin  width  of  A/  =  A/j  as  shown  in  Figure  2.6.  Table  2.2  includes  the  interpolation 
factors  used  to  obtain  the  uniform  bin  width  for  all  subbands.  Note  that  Ni  does  not 
change  as  5  ^  is  not  modified.  "V^fith  the  change  in  the  bin  width,  the  corresponding  data 
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points  in  subbands  S2  -S^  can  be  expressed  as 


(2.7) 

where  is  the  new  number  of  data  points,  A/„  is  the  original  bin  width  of  subband 
signal  ,  for  n  =2, 3 . 6,  and  A/^  is  the  frequency  bin  width  of  subband  signal  5^ . 


Ni 

N4 

^5 

A/  =  A/i 


J 


0  0.51  2.192  40.55  121.6  156.9  243.1  1000  MHz 

Figure  2.6  Conversion  to  Uniform  Bin  Width  by  Interpolation 


Table  2.2:  Subband  Frequency  Ranges  and  Interpolation  Factors  of  Test  Data 


Subband 

Signal 

Frequency  Band 
(Hz) 

Interpolation 

Factor 

Si 

0.510e6  -  2.192e6 

1 

S2 

2.192e6  -  40.55e6 

43/10 

S3 

40.55e6-  121.6e6 

80 

'S4 

121.6e6-159.6e6 

224 

S5 

159.6e6  -  243.le6 

308 

Se 

243.1e6-1000e6 

478 

The  total  data  number  of  points  in  the  frequency  domain,  ,  after  the  bin  width 
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conversion  in  the  subbands,  is  given  by 


^total 


=  Ni  +  N^+ 


(2.8) 


The  bin  width  in  the  frequency  domain  is  related  to  the  sampling  frequency,  f^,  and  the 
total  number  of  points  as  follows; 


A/ 


2A^total 


1 

^A'totalA^ 


(2.9) 


where  At  is  the  corresponding  sampling  period  in  the  time  domain.  The  sampling  period 
or  the  sampling  frequency  can,  therefore,  be  determined  as 


At  = 


1 

2A^totalA/ 


or  f,  =  2A^totalA/ 


(2.10) 


D.  TIME-DOMAIN  REPRESENTATION  OF  THE  CSW  WAVEFORM 


The  time-domain  representation  of  the  CSW  data  can  be  achieved  by  following  the 
steps  outlined  in  the  previous  two  sections.  Given  the  CSW  data  measured  in  six  subbands 
at  different  bin  widths,  we  first  perform  the  bin  width  conversion  to  obtain  a  uniform  bin 
width  representation  as  shown  in  Figure  2.6.  From  Equation  2.10,  we  have  the 
corresponding  samplingperiod. 

To  transform  the  frequency  domain  CSW  data  into  the  time  domain  form,  the 
approach  discussed  in  section  B  is  used.  Following  on  the  lines  of  Figure  2.3,  the  two- 
sided,  zero-padded  CSW  spectrum  is  obtained  as  shown  in  Figure  2.7.  An  inverse  Fourier 
transform  is  then  performed  on  the  spectral  data  to  obtain  the  corresponding  time-domain 
signal. 
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Figure  2.7  Time-Domain  Representation  of  the  CSW  Waveform:  (a)  Mirror  Image  CSW 

Spectra  and  (b)  Zero-Padded  CSW  Spectra. 


Two  methods  were  used  to  produce  the  time-domain  representation  of  the  CSW 

waveform:  Frequency  Domain  interpolation  Method  and  Time  Domain  interpolation 

Method.  The  frequency  domain  interpolation  method  produces  a  time  domain  CSW 

waveform  in  the  following  manner  (see  Figure  2.8): 

Step  1.  Perform  bin  width  conversion  on  the  entire  frequency  domain  data,  so  the 
frequency  domain  samples  of  all  subbands  are  equally  spaced. 

Step  2.  Take  the  inverse  Fourier  transform  of  the  entire  frequency  domain  CSW 
waveform. 
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Figure  2.8  Schematic  Diagram  to  Realize  the  CSW  in  the  Time  Domain  by  First 
Interpolating  the  Subband  Data  in  the  Frequency  Domain. 


The  time  domain  interpolation  method  produces  a  time  domain  CSW  waveform  in  the 
following  maimer  (see  Figure  2.9): 

Step  1.  Take  the  inverse  Fourier  transform  of  each  subband  without  changing  the 
bin  width. 

Step  2.  Interpolate  subband  CSW  waveforms  in  the  time  domain. 

Step  3.  Sum  together  to  obtain  the  time  domain  CSW  waveform. 

The  interpolation  of  the  subband  signals  in  the  time  domain  is  performed  such  that 
Equation  2. 10  is  satisfied.  The  interpolation  factors  used  in  the  time  domain  are  same  as 
those  in  the  frequency  domain  (see  Table  2.2). 
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Figure  2.9  -Schematic  Diagram  to  Realize  the  CSW  in  the  Time  Domain  by  Interpolating 

the  Subband  Data  in  the  Time  Domain. 

Figures  2.10  and  2.11  show  the  plots  of  the  time  domain  CSW  waveform  obtained 
from  the  two  methods.  In  the  frequency  domain,  all  components  below  80  MHz  have  been 
discarded,  which  mainly  leaves  four  subbands:  part  of  S3  and  S4  -Sg  (see  Figure  2.5). 

Figure  2. 12(a)  shows  the  magnitude  spectrum  of  the  original  CSW.  Figure  2. 12(b) 
shows  the  CSW  representation  in  the  time  domain  by  interpolating  in  the  frequency 
domain  first  and  then  transforming  it  to  the  time  domain.  The  result  of  inteipolating  the 
subband  signals  in  the  time  domain  is  shown  in  Figure  2.12(c).  The  method  of 
interpolation  of  frequency  samples  exhibits  some  spurious  time  domain  values  at  the  tail 
end  of  the  signal;  however,  the  spectral  representation  of  this  method  is  slightiy  better  than 
the  time  domain  sample  interpolation.  In  the  rest  of  the  thesis,  the  time  domain  sample 
interpolation  is  used  because  the  waveform  so  obtained  is  cleaner  than  the  other. 
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Figure  2.10  Frequency  Domain  Interpolation:  (a)  83:86  -  121.6MHz.  (b)  84:121.6  - 
156.9MHz.  (c)  85:  156.9  -  243.1MHz.  (d)  85:  243.1  -  lOOOMHz.  (e)  8ynthesis  8ignal  of 

(a)-(d). 
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Figure2.11  Time  Domain  Interpolation:  (a)  83:80-  121.6MHz.  (b)  84:121.6-  156.9MHz. 
(c)  85:  156.9  -  243.1MHz.  (d)  85:  243.1  -  lOOOMHz.  (e)  8ynthesis  8ignal  of  (a)  -  (d). 
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time  (sec) 


frequency  (Hz) 


Figure  2.12  Comparison;  (a)  Original  CS W  waveform  (Recorded  in  Frequency  Domain) 
(b)  Frequency  Domain  Interpolation  Method  (c)  Time  Domain  Interpolation  Method  (All 
Components  below  80  MHz  Are  Discarded). 


E.  SYNTHESIS  OF  THE  WIDEBAND  WAVEFORM  BY  CONCATENATION 


In  this  section,  we  present  methods  to  synthesize  a  single  wideband  signal  from 
HPW  and  CSW  signals  originating  from  the  same  test  point  with  overlapped  frequency 
spectra  in  the  region  80-100  MHz.  Figure  2.13  shows  a  typical  HPW  waveform 
(measured  in  the  time  domain)  and  its  magnitude  spectrum. 


(b)  frequency  (Hz) 


Figure  2.13  HPW  Waveform  (TAIOlOOl):  (a)  Time  Domain  Representation  (b) 

Magnitude  Spectrum. 


Due  to  the  dynamic  range  and  analog  bandwidth  limitation  of  the  HPW  data,  the 
frequency  content  of  the  HPW  waveform  is  considered  reliable  only  up  to  100  MHz.  Thus, 
even  though  higher  frequency  data  are  available,  our  frequency  range  of  interest  is  limited 
to  from  about  0.5  MHz  to  100  MHz.  The  HPW  data  were  measured  using  a  sampling  rate 
of  /j  =  2.1394  X  10^ .  The  bin  width  of  the  HPW  waveform  is  Afj  =  5.223  X  10^  Hz. 
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Our  goal  is  to  merge  the  low  frequency  HPW  and  the  high  frequency  CS  W  to  obtain  a  single 
waveform  covering  the  frequency  range  from  0.5  MHz  to  1  GHz.  The  CSW  phase 
measurements  may  be  inaccurate  below  80  MHz.  This  prevents  our  obtaining  the  low 
frequency  waveform  from  CSW  data  alone.  Figure  2.14  shows  the  magnitude  spectra  of 
both  HPW  (TAIOlOOl)  and  CSW  (LAIOlOOl)  waveforms  which  were  collected  from  the 
same  test  point  (lOlOOl)  in  parallel  (A)  orientation. 


(a) 


frequency  (Hz) 


Figure  2.14  Magnitude  Spectrum  of  (a)  HPW  Waveform  (TAIOlOOl)  and  (b)  CSW 

Waveform  (LAIOlOOl). 


The  HPW  and  CSW  waveforms  are  truncated,  leaving  an  overlapped  region  80  - 
100  MHz.  Figure  2.15  shows  the  truncated  and  overlapped  spectra  of  HPW  and  CSW 
waveforms.  Frequency  components  from  100  to  1000  MHz  of  TAIOlOOl  and  frequency 
components  from  0.51  to  80  MHz  of  LAIOlOOl  were  discarded.  Prior  to  concatenation 
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and  inverse-Fomier  transform  of  the  combined  spectrum,  the  bin  width  of  all  the  subbands 
including  that  of  the  HPW  have  to  be  the  same.  Table  2.3  shows  the  frequency  ranges,  bin 
width  and  interpolation/decimation  factors  needed  to  achieve  a  imiform  bin  width.  Then 
a  two-sided  spectrum  is  formed,  and  the  time  domain  representation  of  the  combined  signal 
is  obtained  by  the  inverse-Fourier  transform  of  this  spectrum. 
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Figure  2.15  Subband  Spectra:  (a)  HPW  Waveform  (TAIOlOOl)  (b)  CSW  Waveform 

(LAIOlOOl). 


Prior  to  combining  the  CSW  and  HPW  waveforms,  the  CSW  data  is  normalized  by 
multiplying  by  a  constant: 


(2.11) 
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Table  2.3:  Subband  Frequency  Ranges,  Bin  Width  and  Interpolation/Decimation 

Factors  of  Test  Data 


Subband 

Signal 

Frequency  Range 
(Hz) 

Bin  Width 
(Hz) 

Bin  Width  Conversion 
(Up/Down)  Factor 

0-  100e6 

5.223 le5 

1/1 

Si 

0.510e6  -  2.192e6 

0.0847e5 

1/62 

S2 

2.192e6  -  40.55e6 

0.365  le5 

3/43 

S3 

40.55e6  -  121.6e6 

6.7375e5 

9/7 

S4 

121.6e6-159.6e6 

19.500e5 

15/4 

S5 

159.6e6-243.le6 

26.060e5 

5/1 

Se 

243.1e6-1000e6 

40.476e5 

31/4 

J[  VQOMHz  j  lOOMHz 

52  Tiif),  ^  ^2^(/)»  and  T(f)  represents  the 

f=  80MHz  Hf  =  80MHz 

HPW  waveform  and  L(/)  the  CSW  waveform. 

The  combined  waveform,  which  covers  the  entire  frequency  band  of  the  given 
HPW  and  CSW  waveforms  together,  is  obtained  in  two  different  ways  as  shown  Figure 
2. 16.  In  the  first  method,  the  HPW  spectrum  up  to  80  MHz  is  concatenated  with  the  CSW 
spectrum  (80  MHz  to  1  GHz).  In  the  second  method,  the  HPW  spectrum  up  to  100  MHz 
is  concatenated  with  the  CSW  spectrum  (100  MHz  to  1  GHz).  Thus,  the  overlapped 
region  is  covered  by  the  components  of  either  HPW  or  CSW.  Later,  another  approach  to 
combine  the  spectral  components  over  the  overlapped  region  is  detailed. 

Before  the  spectral  concatenation  is  carried  out,  the  HPW  and  CSW  waveforms 
have  to  be  time  aligned.  From  Figures  2.11(e)  and  2.13(a),  we  notice  that  the  HPW 
waveform  starts  to  build  up  around  0.2  p r  while  the  CSW  begins  aroimd  the  origin.  Since 
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the  HPW  is  measured  in  the  time  domain,  the  HPW  waveform  is  used  as  a  reference,  and 
the  CSW  waveform  is  time  shifted  to  align  with  the  HPW  waveform. 


t 
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(a)  Concatenation  at  80  MH& 


1000  MHz 
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1000  MHz 


Figure  2.16  Concatenation  of  HPW  and  CSW  spectra;  (a)  80  MHz  and  (b)  100  MHz. 

The  results  of  the  two  methods  are  shown  in  Figure  2.17  and  2.18.  Both  plots  look 
very  similar.  The  sharp  tnmcation  of  spectra  as  described  in  the  above  may  lead  to  some 
spurious  components  in  the  time  domain  unless  the  amplitude  of  the  spectrum  is  small 
enough  at  the  truncation  frequency  [Ref.  6].  In  Figure  2.17  and  2.18,  however,  no  such 
components  are  visible.  Nevertheless,  any  effects  of  truncation  can  be  minimized  by 
making  a  gentler  transition  to  zero  at  the  cutoff  frequency  through  a  proper  filter  function. 
A  possible  disadvantage  of  the  use  of  the  filter  fimction  is  that  there  is  an  overlap  at  the 
cutoff  frequency  which  should  be  investigated.  The  use  of  a  filter  bank,  instead  of  the 
frequency  truncation  as  above,  will  be  addressed  in  Chapter  IH. 
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10®  10®  10^  10®  10® 


(b)  frequency  (Hz) 

Figure  2.18  Concatenation  at  100  MHz:  (a)  Time  and  (b)  Frequency  Representation. 
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III.  MULTIRATE  FILTER  BANKS 


The  area  of  multirate  digital  signal  processing  is  an  important  part  of  the  modem 
(digital)  communications  theory  in  which  digital  transmission  systems  are  required  to 
handle  data  at  several  sampling  rates  [Ref.  7].  This  technique  finds  applications  in  speech 
and  image  compression,  digital  audio,  statistical  and  adaptive  signal  processing,  and  in 
many  other  fields.  It  also  can  be  applied  to  certain  classes  of  time-frequency 
representations  such  as  the  short-time  Fourier  transform  and  the  wavelet  transform,  which 
are  useful  in  analyzing  the  time- varying  nature  of  wideband  transient  signal  spectra  [Ref. 
3]. 

Multirate  filter  banks  provide  flexibility  in  the  analysis  of  signals  by  dividing  the 
original  signal  into  smaller,  more  manageable  subbands.  According  to  the  multirate  filter 
bank  theory  [Ref.  8],  there  is  no  need  to  satisfy  the  sampling  theorem  on  a  channel  by 
channel  basis  because  it  is  sufficient  to  satisfy  it  on  the  sum  of  the  channels.  Therefore, 
the  impractical  ideal  bandpass  filters  are  not  necessary  anymore,  and  a  new  theory  can  be 
developed  which  looks  at  the  channels  simultaneously  instead  of  regarding  each  one 
separately  as  in  single  filter  signal  processing. 

A.  THEORETICAL  BACKGROUND 

1.  Perfect  Reconstruction  (PR)  Filter  Bank 

It  is  well  known  that  the  quadrature  mirror  filter  (QMF)  bank  can  be  used  in  both 
analysis  and  synthesis  banks  to  produce  a  perfect  reconstruction  synthesis  signal  [Ref.  3]. 
The  most  weU-known  type  of  multirate  filter  bank  is  the  M-chaimel  maximally  decimated 
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QMF  bank,  which  is  shown  in  Figure  3.1. 


The  reconstructed  signal,  x{n) ,  in  general  suffers  from  aliasing  error  because  the 
analysis  bank  filters,  H^(z) ,  that  precede  the  decimators  are  not  ideal.  However,  in 
practice,  for  a  given  set  of  analysis  filters,  ,  the  synthesis  filters,  Fj^iz) ,  can  be 
chosen  to  reduce  the  effect  of  the  aliasing  caused  by  the  decimation  operation. 

The  most  general  expression  for  X(z)  is  of  the  form  [Ref.  3] 

M-l  M-l 

=  5  E  E  HJ^zW-‘)Ft(z)  (3.1) 

/  =  0  jt  =  0 

where  X(zW  =  X(ze^  jnl)/M^ ,  /  ^  0  represents  the  aliasing  terms.  If  we  consider  a 
two  channel  (M -2)  filter  bank,  this  ejqsression  can  be  rewritten  as 
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(3.2) 


Mz)  =  \[HQiz)FQiz)  +  H^iz)F^(z)]Xiz) 

+  ^[HQ(-z)FQiz)  +  H^(-z)Fi(z)]Xi-z) 

where  X(-z)  corresponds  to  on  the  unit  circle.  The  frequency  response  of 

the  two  channel  filter  bank  is  represented  by  the  first  term  in  Equation  3.2  while  the  second 
term  represents  the  aliasing.  To  obtain  aliasing  cancellation,  the  analysis  and  systhesis 
filters  must  be  selected  so  that  the  quantity  of  the  second  term  is  zero: 

Fo(z)  =  H^i-z) 

r  r  \  u  /  \ 

Fi(z)  =  -Hq{-z). 

Therefore,  once  the  analysis  filters,  Hq{z)  and  H^iz),  are  given,  it  is  possible  to 
completely  cancel  aliasing  by  the  choice  of  matching  systhesis  filters,  Fq(z)  and  Fj(z) , 
which  are  lowpass  and  highpass  respectively.  Using  Equation  3.3  for  the  aliasing 
cancellation  conditions.  Equation  3.2  becomes 

X(z)  =  i[/fo(z)Fo(z)  +  //i(z)Fi(z)]Z(z)  =  nz)X{z)  (3.4) 

1 

where  r(z)  =  -  [//q(z)Fq(z)  +  (z)F2(z)] .  If  r(z)  satisfies  the  conditions  of  a  stable 

all-pass  and  a  linear-phase  FIR  function,  then  this  system  is  free  from  amplitude  distortion 
and  phase  distortion,  respectively.  Such  QMF  banks  which  are  free  from  aliasings 
amplitude  distortion  and  phase  distortion  are  called  perfect  reconstruction  (PR)  banks  and 
satisfy 


X(z)  =  cz  '^X(z)  or  x(n)  =  cx(n  -  Mq)  (3.5) 

where  c  is  an  arbitrary  nonzero  constant  and  /Iq  is  a  positive  integer  [Refs.  3  and  5]. 
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2.  Perfect  Reconstruction  Non-uniform  Filter  Bank  With  Rational  Sampling 


Factors 


Figure  3.2  shows  a  general  non-uniform  filter  bank  with  rational  decimation  factors. 
In  this  system,  an  input  signal  can  be  arbitrarily  decomposed  into  a  number  of  different  non- 
uniform  frequency  subbands.  For  a  perfect  reconstruction  with  maximally  decimated 
system,  the  decimation  and  interpolation  factors,  Pj^  and  Q^.,  k  =  0,  ...,M  -  1 ,  satisfy  the 
following  relationship: 


ifc  =  0 


(3.6) 


Qm- 


M~1 


M-1 


KKIO 


Analysis  Filters 


_  Qm-1 

-C>00 


x(n) 


Synthesis  Filters 


Figure  3.2  An  Af -band  Analysis-Synthesis  Non-Uniform  Filter  Bank  with  Rational 

Decimation  Factors.  After  [Ref.  9]. 
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B.  FILTER  BANK  DESIGN 


1.  FIR  Power  Symmetric  Filter  Bank 


Using  Equation  3.3,  the  aliasing  cancellation  condition.  Equation  3.4  can  be 
rewritten  as 


1(2)  =  i[Ho(2)//i(-z)  -  (3.7) 

Assume  that  Hq{z)  is  power  symmetric,  that  is 

Hq{z)H^{z)  +  Hq{-z)Hq{-z)  =  \  (3.8) 


where  the  represents  “paraconjugation”  (defined  by  Vaidyanathan  [Ref.  3]  as 
conjugation  of  the  function’s  coefficients  and  replacement  of  z  with  z~^ ).  With  Hq{z) 
satisfying  Equation  3.8,  the  filter,  //^(z) ,  is  chosen  as 

t 

H^{z)  =  -z-^Hq{-z)  (3.9) 


where  N  is  some  odd  number.  The  remaining  filters  can  be  derived  as  follows: 


Fo(z)  =  z-%(z) 
Fi(z)  =  z~^Hi{z) 


(3.10) 


Now  only  the  filter,  Hq{z)  ,  remains  to  be  designed.  The  power  symmetric  property  means 
that  the  zero-phase  filter,  H{z)  =  Hq{z)Hq{z)  ,  is  a  half-band  filter.  The  design  steps  for 
the  real  coefficient  case,  where  h^in)  is  real,  are  as  follows: 

Step  1.  Design  a  zero-phase  FIR  half-band  filter,  G(z)  =  _  ^8(n)z~" . 
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Step  2.  Define  H{z)  =  G{z)  +  6 ,  where  6  is  the  peak  stopband  ripple  of  G(e^^) 
to  ensure  that  is  nonnegative. 

Step  3.  Compute  Hq(z)  as  a  spectral  factor  of  H(z)  [Ref  3]. 

Figure  3.3  shows  an  example  of  perfect  reconstruction  filter  bank  design.  Since  the 
zero-phase  FIR  half-band  filter,  G(e^^) ,  is  real,  is  obtained  by  lifting  the  response 

G(e^^)  by  6 .  This  is  demonstrated  in  Figure  3.3(a). 

The  zero-phase  FIR  filter,  G(z) ,  is  designed  using  MATLAB  function  “REMEZ” , 
which  uses  the  McClellan-Parks  algorithm.  The  noimegative  zero-phase  filter,  Hiz) ,  is 
also  demonstrated  in  Figure  3.3(a).  It  is  clear  that  ^0  for  all  co ,  so  it  is  possible 

to  find  a  spectral  factor,  Hq{z)  ,  of  Hiz) . 

The  frequency  response  and  the  zero  location  of  perfect  reconstruction  filter  banks, 
HQiz)  and  Hiiz)  (filter  length  of  9,  M=2),  are  shown  in  Figures  3.3(b)  and  3.3(c), 
respectively.  The  passband  gain  is  unity,  the  stopband  gain  is  8.5dB,  and  the  cutoff 
frequen<^  is  Ji/2  radians.  As  demonstrated  in  Figure  3.3(c),  the  spectral  factor,  HQiz) , 
of  the  filter,  Hiz) ,  is  not  unique  because  the  zeros  of  Hiz)  can  be  grouped  into  those  of 
HQiz)  and  HQiz)  in  many  ways.  Table  3.1  shows  the  filter  coefficients  of  Ao(«)  and 
hiin). 
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Normalized  Frequency 


o:  zeros  of  Ho(z),  x:  zeros  of  Hi(z) 


(b)  (c) 

Figure  3.3  Perfect  Reconstruction  FIR  Power  Symmetric  Filter  Bank  (Filter  C)rder=9, 
M=2):  (a)  Zero-Phase  FIR  Half-band  Filter  H(z),  (b)  Perfect  Reconstruction  FIR  Filter 
Bank  and  (c)  Zeros  of  Analysis  Filter  Ho(z)  and  Hi(z). 
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Table  3.1:  Perfect  Reconstruction  Filter  Coefficients  (Filter  Order  =  9) 


n 

ho(n) 

hi(n) 

0 

0.5000 

0.5000 

1 

0.5323 

-0.5323 

2 

0.2495 

0.2495 

3 

-0.1059 

0.1059 

4 

-0.1417 

-0.1417 

5 

0.0422 

-0.0422 

6 

0.1001 

0.1001 

7 

-0.0152 

0.0152 

8 

-0.1202 

-0.1202 

9 

0.1129 

-0.1129 

2.  Perfect  Reconstruction  Uniform  and  Non-Uniform  Filter  Banks 

Both  uniform  and  non-uniform  filter  banks  are  designed  in  this  thesis  for  the 
analysis  and  synthesis  filter  banks  in  order  to  obtain  the  desired  subband  signals.  The 
proper  combination  of  these  filter  banks  provides  flexibility  in  generating  subband  signals. 
The  desired  portions  of  a  given  signal  can  be  filtered  into  narrow,  non-uniform  subbands, 
and  the  remaining  portions  of  the  si^al  can  be  filtered  into  larger,  uniform  subbands  by 
selecting  proper  non-uniform  filter  banks  and  uniform  filter  banks  respectively. 

In  the  case  of  the  uniform  filter  banks,  all  of  the  subbands  of  each  filter  have  the 
same  bandwidth.  A  perfect  reconstruction,  maximally  decimated  analysis-synthesis 
system  based  on  imiform  filter  banks  is  designed  as  shown  in  Figure  3.4.  The 
corresponding  synthesis  filter  bank  can  easily  be  designed  from  the  analysis  filter  bank. 
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Figure  3.4  Perfect  Reconstruction  Uniform  Analysis  Filter  Bank  (Filter  Order=35,  M=5); 
(a)  Block  Diagram  and  (b)  Frequency  Response. 


The  frequency  response  of  a  non-uniform  filter  bank  with  integer  downsampling 
rate  is  shown  in  Figure  3.5.  This  non-uniform  system  can  have  a  munber  of  different  forms 
and  structures.  However,  for  non-uniform  filter  banks,  the  decimation  factors,  Dj^, 
k  =  0, . . .  ,M  -  1 ,  have  the  following  constraint 


M-l 


=  1 


jt  =  0 


(3.11) 


The  synthesis  filter  banks  of  both  uniform  and  the  non-uniform  kind  also  can  easily  be 
obtained  from  Equation  3.9  in  that 

F^(z)  =  (3.12) 
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where  F^(z)  and  Hj^iz)  are  the  synthesis  and  analysis  filter  bank,  respectively, 

^  =  0,  1 . 


(a)  (b) 

Figure  3.5  Perfect  Reconstruction  Non-Uniform  Analysis  Filter  Bank  (Filter  Order=35, 
M=6):  (a)  Block  Diagram  and  (b)  Frequency  Response. 

3.  Perfect  Reconstruction  Non-uniform  Filter  Bank  With  Rational  Sampling 
Factors  (Direct  Method) 

A  direct  method  is  presented  in  this  thesis  to  design  the  non-uniform  filter  banks 
with  rational  decimation  factors.  Especially,  a  two-band  non-uniform  filter  bank  with 
rational  decimation  factors  of  5/4  and  5/1,  so  called  (4/5,  1/5)  system,  is  considered  to 
discuss  the  design  process  to  obtain  the  desired  spectral  detail  as  shown  in  Figure  3.6. 
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Figure  3.6(b)  shows  the  frequency  response  of  the  analysis  filter  banks.  According  to  non- 
uniform  filter  bank  theory  [Ref.  9],  a  (4/5, 1/5)  system  can  also  be  designed  by  combining 
the  first  four  bands  of  a  uniform  five-band  system  (see  Figure  3.4). 


x(n) 


<2oU 

x(n) 


(a) 


Figure  3.6  Filter  Bank  with  Rational  Sampling  Factors  (Filter  Order  =  35):  (a)  A  block 
Diagram  and  (b)  The  Frequency  Response  of  (4/5, 1/5)  System. 
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4.  IVee  Structured  Filter  Bank 


Both  uniform  and  nun-uniform  filter  banks  were  arranged  into  a  tree  structured 
analysis  and  synthesis  filter  banks  to  synthesize  two  signals  obtained  from  different 
measurement  devices.  In  this  scheme,  each  input  signal  is  divided  into  smaller,  more 
manageable  subband  signals  which  retain  the  frequency  characteristics  of  the  original 
signal.  They  are  then  summed  to  obtain  a  wideband  signal,  which  covers  the  entire 
frequency  area  and  retains  the  basic  spectral  properties  of  the  original  signals. 

The  overlapped  portions  of  the  two  signals  are  not  handled  separately  in  the  section. 
The  study  of  overlapped  portion  will  be  presented  in  Chapter  IV.  As  mentioned  in  Chapter 
n,  frequency  concatenation  methods  were  applied  in  this  section  using  tree  structured  filter 
banks  to  produce  a  desired  signal.  Equation  2.11  is  applied  to  weigh  the  high  frequency 
portion  of  CSW  data. 

After  decimation  in  a  subband,  a  minimum  data  length  of  three  times  the  filter  order 
is  required  in  order  for  the  filtoing  operation  to  properly  occur  at  the  next  tree  level;  a 
proper  data  length  is  required  to  minimize  the  end-point  transients  that  occur  during  the 
filtering  process  [Ref.  2].  In  the  overlapping  frequency  band  of  80  -  100  MHz,  for 
example,  the  input  is  decimated  by  a  factor  of  50,  which  means  that  only  82  samples  are 
present  in  this  subband.  For  a  filter  order  of  35,  at  least  105  points  are  required  in  this 
subband.  To  overcome  this  problem,  a  “zero-padding”  method  [Ref.  2]  was  implemented 
whenever  the  data  points  fell  short  of  the  required  length. 

Four  tree-structured  filter  banks  were  designed  to  obtain  a  synthesis  signal; 
concatenation  at  80  MHz,  concatenation  at  100  MHz,  concatenation  at  80  MHz  with  (4/5, 
1/5)  system,  and  symmetric  overlap  at  90  MHz.  The  perfect  reconstruction  FIR  power 
symmetric  filter  bank  was  applied  in  the  case  of  overlap  at  90  MHz.  The  FIR  power 
symmetric  filter  bank  provides  perfect  reconstruction  of  the  data,  so  it  performs  better  than 
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truncation  by  providing  a  smooth  transition  between  the  HPW  and  the  CSW  spectra. 

Figures  3.7-3.10  show  the  schematic  diagrams  for  the  four  filter  banks.  The  first 
filter  bank  takes  HPW  data  up  to  80  MHz  and  CSW  data  from  80  MHz  to  1000  MHz,  and 
the  second  filter  bank  takes  HPW  data  up  to  100  MHz  and  CSW  data  from  100  MHz  to 
1000  MHz  to  obtain  a  wideband  transient  waveform.  The  tree  structured  filter  bank  using 
rational  decimation  factor  of  5/4  and  5/1  is  shown  in  Figure  3.9.  The  fourth  filter  bank 
takes  HPW  data  up  to  90  MHz  and  CSW  data  from  90  MHz  to  1000  MHz  as  shown  in 
Figures  3.10  and  3.11;  the  frequency  response  of  the  power  symmetric  filter  bank  around 
90  MHz  is  shown  in  Figure  3.12. 
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Figure  3.7  Analysis  and  Synthesis  Tree  Structured  Filter  Bank  with  Concatenation  at  80 

MHz. 
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Figure  3.8  Analysis  and  Synthesis  Tree  Structured  Filter  Bank  with  Concatenation  at  100 

MHz. 
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Figure  3.9  Analysis  and  Synthesis  Tree  Structured  Filter  Bank:  Concatenation  at  80MHz 
with  Rational  Decimation  Factor  of  5/4  and  5/1. 


40 


Figure  3.10  Analysis  and  Synthesis  Tree  Structured  Filter  Bank  with  Power  Symmetric 


Overlap  at  90  MHz. 
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frequency  (MHz) 

Figure  3.12  Frequency  Response  of  Tree-Siructured  Filter  Bank  with  Symmetric  Overlap 

at  90  MHz  (Filter  Order  =  35). 
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Figure  3.13  shows  the  results  obtained  from  the  tree  structured  filter  bank  with 
concatenation  at  80  MHz  for  test  point  lOlOOl.  Figure  3.14  shows  the  results  obtained 
from  the  tree  structured  filter  bank  with  concatenation  at  100  MHz  for  the  same  test  point. 
The  results  obtained  from  tree  structured  filter  bank  (concatenation  at  80  MHz)  with 
rational  decimation  factor  of  5/4  and  5/1  are  shown  in  Figure  3.15.  Figure  3.16  shows  the 
results  obtained  from  the  tree  structured  filter  bank  with  symmetric  overlap  at  90  MHz. 

The  plots  in  Figures  3.13-3.16  look  very  similar  expect  for  slight  differences  at  the 
concatenation  point. 


Waveform  and  (b)  Magnitude  Spectrum. 
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(b)  frequency  (Hz) 


Figure  3,14  Tree-Structured  Filter  Bank  with  Concatenation  at  100  MHz:  (a)  Time 
Domain  Waveform  and  (b)  Magnitude  Spectrum. 
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(b)  frequency  (Hz) 

Figure  3.15  Tree  Structured  Filter  Bank  (Concatenation  at  80  MHz)  with  (4/5  and  1/5) 
System:  (a)  Time  Domain  Waveform  and  (b)  Magnitude  Spectrum. 
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Figure  3.16  Tree-Structured  Filter  Bank  with  Symmetric  Overlap  at  90  MHz:  (a)  Time 
Domain  Waveform  and  (b)  Magnitude  Spectrum. 
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rV.  WIDEBAND  ANALYSIS  RESULTS 


The  results  of  systhesizing  the  waveform  data  using  bandpass  IFFT  and  tree- 
structured  filter  bank  are  presented  in  this  chapter.  In  Chapters  II  and  IH,  the  bandpass 
IFET  filter  bank  method  and  tree-structured  filter  bank  method  with  firequency 
concatenation  scheme  were  applied  to  combine  HPW  and  CSW  waveforms,  respectively. 
In  this  chapter,  the  sinusoidal  modeling  algorithm,  the  bandpass  IFFT,  and  the  tree- 
structured  filter  bank  techniques  will  be  used  to  combine  four  different  waveforms  at  a 
single  test  point  to  produce  a  synthesized  test  point  signal. 

A.  SINUSOIDAL  MODELING  ALGORITHM  (SMA) 


The  main  idea  of  this  technique  is  that  any  waveform  can  be  represented  as  a  sum 
of  sine  waves.  The  sinusoidal  modeling  technique  implemented  in  this  chapter  was 
originally  developed  by  McAxilay  and  Quatieri  [Ref.  10]  for  speech  signal  processing, 
which  is  transient  in  nature  and  stationary  for  only  a  short  interval  of  time. 

The  block  diagram  illustrating  the  SMA  is  shown  in  Figure  4.1.  In  the  analysis 
section,  the  input  data  are  sectioned  off  into  equal  length  frames  and  weighted  with 
Hamming  window.  Frames  are  formed  allowing  for  overlap  of  data.  The  short-time 
Fourier  transform  (STFT)  approach  is  used  to  obtain  the  amplitudes,  frequencies  and 
phases  of  the  sine  waveforms  in  each  spectral  bin. 

The  frequency  values  are  determined  by  first  locating  the  peaks  in  the  periodogram, 
and  then  matching  these  with  peaks  in  adjacent  frames.  These  peaks  need  not  occur  at  the 
same  spectral  bin  but  are  constrained  to  lie  within  some  matching  interval.  In  the  synthesis 
part  (see  Figure  4.1(b)),  a  corresponding  linear  amplitude  interpolation  and  a  phase 
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unwrapping  and  interpolation  between  matching  peaks  are  performed.  For  each  set  of 
matching  peaks,  a  sine  wave  is  generated  and  amplitude  modulated.  The  resulting  sine 
waves  are  then  summed  to  produce  the  synthesized  waveform.  For  more  details  on  this 
algorithm,  see  McAulay  and  Quatieri  [Ref.  10]. 


(b) 


Figure  4.1  Block  Diagram  of  the  Sinusoidal  Modeling  Algorithm:  (a)  Analysis  and  (b) 

Synthesis  System. 
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In  this  thesis,  we  have  used  the  following  parameters  for  the  sinusoidal  modeling 
algorithm:  a  FFT  size  of  512  and  a  frame  length  of  256  were  chosen  to  obtain  proper  bin 
width  and  frequency  resolution;  an  overlap  interval  of  128  was  chosen  for  50%  overlap;  a 
matching  window  size  of  3  was  chosen  to  ensure  smooth  frequency  track;  and  a  peak 
threshold  of  30  dB  was  chosen  to  reject  weak  spectral  components  and  others  due  to  noise 
[Ref.  1]. 

B.  COMBINATION  OF  WAVEFORMS  FROM  DIFFERENT  ORIENTATIONS 

1.  Data  Characteristics 

In  Chapters  n  and  IH,  to  test  the  various  algorithms,  only  the  data  from  A  (parallel) 
orientation  were  considered.  The  B  (perpendicular)  orientation  data  will  now  be  analyzed 
to  combine  them  into  one  waveform.  There  are  two  different  sets  of  waveform  data  at  each 
test  point:  HPW  parallel  and  HPW  perpendicular;  and  CSW  parallel  and  CSW 
perpendicular. 

Figures  4.2  and  4.3  show  the  four  waveforms  measured  at  test  point  lOlOOl .  Figure 
4.2  shows  the  two  HPW  waveform  and  their  spectra.  The  B  (perpendicular)  orientation 
spectrum  shows  a  large  low  frequency  coupling  below  the  6  MHz  frequency  area.  The 
frequency  content  in  the  frequericy  band  of  6  - 100  MHz  is  different  for  the  two  waveforms. 
The  time  resolution  and  the  sampling  frequency  for  both  TAIOlOOl  (parallel)  and 
TBIOKX)!  (perpendicular)  orientation  waveforms  are  given  in  Table  4.1.  The  sampling 
periods  are  slightly  different;  however,  in  this  work,  we  have  treated  them  to  be  the  same; 
we  have  used  At  =  0.4674  x  lO^(sec) ,  which  is  from  TAI01(X)1. 
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Table  4.1:  HPW  Waveforms  at  Test  Point  lOlOOl 


TAIOlOOl 

TBIOlOOl 

Number  of  Data  Point 

4096 

4096 

Sampling  Period,  At  (sec) 

0.4674e-9 

0.4629e-9 

Sampling  Frequency,  (Hz) 

2.1394e9 

2.1603e9 

Bin  Width,  A/  (Hz) 

5.2231e6 

5.2742e6 
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(d)  frequency  (Hz) 

Figure  4.2  HPW  Data;  (a)-(b)  Time  Domain  Waveform  and  Magnitude  Spectrum  of 
TAIOlOOl;  (c)-(d)  Time  Domain  Waveform  and  Magnitude  Spectrum  of  TBIOIOOI. 

Figure  4.3  shows  the  originally  recorded  frequency  domain  CSW  data  and  its  time 
domain  representation  at  test  point  lOlOOl .  The  frequency  content  of  the  two  waveforms 
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Figure  4.3  CSW  Data:  (a)-(b)  Magnitude  Spectrum  and  Time  Domain  Waveform  of 
LAIOlOOl;  (c)-(d)  Magnitude  Spectrum  and  Time  Domain  Waveform  of  LBIOIOOI. 


2.  Combination  of  the  HPW  and  CSW  Waveforms  using  Sinusoidal  Modeling 
Algorithm 

The  two  waveforms  are  analyzed  using  the  short-time  Fourier  transform  (STFT) 
approach,  which  yields  a  set  of  frequencies,  magnitudes  and  phases  for  each  waveform. 
The  two  sets  of  analysis  data  are  input  to  the  peak  picking  block  (see  Figure  4.1(a))  of  the 
SMA  algorithm.  For  a  given  sinusoidal  frequency,  the  peak  picking  algorithm  tries  to 
match  the  spectral  peaks  in  the  two  spectra.  The  spectral  peaks  and  the  corresponding 
phases  are  then  input  to  the  synthesis  part  of  the  algorithm  (see  Figure  4.1(b)). 

Figure  4.4  shows  the  plots  of  TAIOlOOl  (parallel)  and  TBIOlOOl  (perpendicular) 
orientation  waveforms  and  the  combined  waveform  obtained  using  the  sinusoidal  modeling 
algorithm.  Figure  4.5  shows  the  plots  of  LAIOlOOl  (parallel)  and  LBIOlOOl 
(perpendicular)  orientation  waveforms  and  the  combined  waveform  obtained  using  the 
sinusoidal  modeling  algorithm.  Before  the  CSW  waveforms  are  combined,  they  were 
time  shifted.  The  time  shift  was  necessary  to  ensure  the  same  sequence  of  windowing  on 
both  HPW  and  CSW  data.  There  are  fewer  samples  in  the  synthesized  waveform  than 
those  in  the  input  waveforms,  which  is  a  result  of  applying  the  SMA  algorithm. 
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magnitude  magnitude  magnitude 


Figure  4.4  Sinusoidal  Modeling  Algorithm  Applied  to  the  HPW  Waveforms:  (a) 
TAIOlOOl  (b)  TBIOlOOl  and  (c)  Combined  Waveform. 


Figure  4.5  shows  the  plots  of  LAIOlOOl  (parallel)  and  LBIOlOOl  (perpendicular) 
orientation  waveforms  and  the  combined  waveform  obtained  using  the  sinusoidal  modeling 
algorithm. 
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Figure  4.5  Sinusoidal  Modeling  Algorithm  Applied  to  the  CSW  Waveforms:  (a) 
LAIOlOOl  (b)  LBIOlOOl  and  (c)  Combined  Waveform. 
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C.  SYNTHESIS  OF  WIDEBAND  TRANSIENT  SIGNALS  FROM  FOUR 


DIFFERENTLY  ORIENTED  DATA  SETS 


In  this  section,  techniques  developed  in  Chapter  n  and  HI  are  applied  to  combine 
the  HPW  and  the  CSW  waveforms  synthesized  in  the  previous  section.  The  result  is  a 
synthesized  test  point  waveform  obtained  as  a  combination  of  measurements  from  two 
different  orientations  and  two  different  spectral  ranges.  Here,  we  use  the  bandpass  IFFT, 
multirate  filter  banks,  and  the  SMA  algorithm  to  obtain  the  synthesized  test  point 
waveform.  Results  of  the  two  techniques  are  presented  for  the  case  where  the 
concatenation  take  place  at  90  MHz  for  the  bandpass  IFFT  technique  and  the  power 
symmetric  overlap  takes  place  over  a  2  MHz  frequency  range  around  90  MHz  for  the  filter 
bank  technique. 

In  the  first  method,  the  bandpass  IFFT  method  with  spectral  concatenation  at  90 
MHz  is  used.  First,  the  HPW  and  the  CSW  waveforms  from  the  same  orientation  are 
combined  into  a  single  synthesized  waveform  using  the  spectral  concatenation  technique 
(developed  in  Chapter  II)  at  90  MHz.  Then,  the  two  synthesized  waveforms  (from  A  and 
B  orientations)  are  combined  using  the  SMA  algorithm  to  obtain  the  final  synthesized  test 
point  waveform.  In  the  second  method,  the  tree-structured  filter  bank  method  (developed 
in  Chapter  n)  with  power  symmetric  overlap  around  90  MHz  is  used  with  the  same 
sequence  of  operations- as  in  the- first  ihethod. 

In  both  methods,  we  could  have  reversed  the  sequence  of  operations  to  obtain  the 
synthesized  test  point  waveform.  The  two  HPW  waveforms  (from  A  and  B  orientations) 
could  be  combined  into  a  single  HPW  waveform,  and  the  two  CSW  waveforms  (from  A 
and  B  orientations)  could  be  combined  into  a  single  CSW  waveform  using  the  SMA 
algorithm.  Then,  these  two  waveforms  could  be  combined  using  either  the  bandpass  IFFT 
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or  the  filter  bank  technique.  The  results  of  this  approach  (SMA  first)  and  the  one  detailed 
above  (SMA  last)  were  similar.  In  this  section,  we  only  present  the  results  using  the  SMA 
last  approach. 

The  results  of  the  two  methods  are  presented  in  Figures  4.6  and  4.7.  In  each  case, 
the  synthesized  test  point  waveform  has  3712  data  points.  This  waveform  contains  384 
fewer  samples  than  the  input  data  sets,  whose  length  is  4096.  The  difference  in  data 
lengths  between  file  input  and  output  data  sets  is  attributed  to  the  data  point  suppression 
characteristic  of  the  sinusoidal  modeling  algorithm  as  mentioned  earlier.  The  time  and 
frequency  domain  representations  of  the  synthesized  test  point  waveform  obtained  from  the 
two  methods  are  included  in  Figure  4.8,  so  the  two  final  waveform  can  be  visually 
compared. 
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Figure  4.6  Bandpass  IFFT  with  Concatenation  at  90  MHz:  (a)  Synthesized  Time  Domain 
Waveform  of  TAIOlOOl  and  TBIOlOOl;  (b)  Synthesized  Time  Domain  Waveform  of 
LAIOlOOl  and  LBIOlOOl;  and  (c)  Synthesized  Test  Point  Waveform. 
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magnitude  magnitude  magnitude 


Figure  4.7  Filter  Bank  with  Symmetric  Overlap  around  90  MHz:  (a)  Synthesized  Time 
Domain  Waveform  of  TAIOlOOl  and  LAIOlOOl;  (b)  Synthesized  Time  Domain  Waveform 
of  TBIOlOOl  and  LBIOlOOl;  and  (c)  Synthesized  Test  Point  Waveform. 
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Figure  4.8  Synthesized  Test  Point  Waveforms:  (a)  Time-domain  and  (b)  Frequency- 
domain  Representations  using  the  Bandpass  IFFT  with  Concatenation  at  90  MHz;  (c) 
Time-domain  and  (d)  Frequency-domain  Representations  using  the  Filter  Bank  with 

Symmetric  Overlap  around  90  MHz. 
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D.  NORM  ATTRIBUTE  MEASURE  RESULTS 


The  true  test  point  waveform  that  takes  the  different  airoraft  orientations  into 
account  and  covers  the  frequency  range  from  0.5  MHz  to  1  GHz  cannot  be  physically 
measured.  No  single  measurement  device  is  currently  available  that  can  measure  a 
wideband  signal  with  its  spectrum  extending  from  0.5  MHz  to  1  GHz.  The  techniques 
developed  in  the  previous  sections  take  the  four  components  waveforms  measured  imder 
two  different  orientations  and  over  two  different  frequency  ranges  and  produce  a 
synthesized  test  point  waveform.  Since  we  do  not  have  a  true  test  point  waveform 
available,  there  is  no  simple  way  of  determining  if  the  synthesized  test  point  waveforms  that 
we  produced  are  indeed  the  desired  waveforms.  One  may  visually  compare  the  different 
component  waveforms  with  the  synthesized  test  point  waveform  to  derive  some  subjective 
observations.  In  this  section,  we  use  the  norm  attribute  measure  [Ref.  1 1]  in  an  attempt  to 
provide  a  means  of  quantitative  comparison. 

A  norm  attribute,  N,of&  time  domain  waveform,  /(r) ,  is  written  as  [Ref.  1 1] 


iV  =  11/(011^0.  (4.1) 

The  norm  attribute  is  a  scalar  having  the  following  properties: 

N  =  0-  if  and  only  if  /(r)  =  0,  for  all  r  (4.2) 

IU/(0lI  =  Ul  11/(011.  A  is  any  scalar  (4.3) 

||/i(0  +  /2(0||  ^  ||/i(0||  + 1|/2(/)||  (triangle  inequality)  (4.4) 
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Table  4.2  includes  the  definitions  for  a  set  of  five  norm  attributes  which  are  believed  to 
capture  the  important  features  of  electromagnetic  transient  waveforms.  For  a  detailed 
discussion  of  the  five  norm  attributes,  see  Thomas  and  Lubell  [Ref.  11]. 

The  first  norm,  ,  is  the  peak  amplitude  of  the  waveform  (in  this  case  magnitude 
of  current) .  The  second  norm,  A^2 » is  the  peak  value  of  the  derivative  of  the  waveform.  The 
third  norm,  N'^,  is  the  peak  impulse  and  measures  the  peak  charge  for  current 
measmements.  The  fourth  norm,  ,  measures  the  total  area  of  the  waveform,  and  the 
fifth  norm,  ,  measures  the  waveform  energy  [Refs.  1  and  11]. 


Table  4.2:  Norm  Attributes,  After  [Ref.  11] 


Norm  Quantity 

Norm 

Attribute 

Feature  of  f(0 

Ni  =  l/«)L„ 

Peak 

Amplitude 

Amplitude 

N2  = 

Peak 

Derivative 

Variation 

Peak 

Impulse 

Content 

=  ^\fix)\dx 

Rectified 

Impulse 

Content 

0  1/2 

N5  =  [j^f  (x)dx] 

Root  Action 
Integral 

Content 

These  norm  attributes  are  used  to  determine  whether  the  synthesized  test  point 
waveform  retains  the  largest  (or  worst  case)  norm  attributes  firom  the  original  component 
waveforms.  The  norm  attributes  of  the  component  waveforms  and  the  synthesized  test 
point  waveform  for  test  point  lOlOOl  are  listed  in  Table  4.3  and  4.4  for  the  bandpass  IFFT 
method  and  the  filter  bank  metiiod,  respectively. 
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From  Table  4.3,  both  methods  produce  the  norm  attributes  that  are  very  close. 
Except  for  A^4 ,  the  other  norm  values  of  the  synthesized  test  point  waveform  are  close  to 
the  maximum  values  of  the  component  waveforms.  Comparison  over  the  overlapped 
region,  from  Table  4.4,  indicates  that  the  filter  bank  technique  produces  norm  values  that 
are  closer  to  the  maximum  norm  values  of  the  component  waveforms. 


Table  4.3;  Norm  Attribute  Comparison  of  the  Component  and  Synth^ized 

Waveforms  (Test  Point  101001) 


Norm 

Attribute 

Cqn^nentWavB&ams 

Syofliesized  Test 
Point  Wav^<«ms 

TA 

01001 

TB 

01001 

LA 

01001 

LB 

01001 

Bandpass 

IFFT 

Filter 

Bank 

Peak 

Amplitude 

0.5760 

e-6 

0.5480 

e-6 

0.3900 

e-6 

0.3890 

e-6 

Peak 

Derivative 

0.2125 

c+8 

0.2139 

6+8 

4.1121 

e+8 

4.9478 

e+8 

2.7869 

6+8 

2.3438 

6+8 

Peak 

Impulse 

0.0900 

e-6 

0.0006 

e-6 

m 

Rectified 

Impulse 

0.6250 

e-6 

0.0780 

e-6 

0.0410 

e-6 

0.0040 

e-6 

Root  Action 
Integral 

0.2837 

0.3382 

0.4460 

0.5682 

0.5022 

0.5032 
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Table  4.4:  Norm  Attribute  Comparison  of  the  Component  and  Synthesized 
Waveforms:  Overlapped  Area  from  80  to  100  MHz  only  (Test  Point  lOlOOl) 


Norm 

Attribute 

Componeaat  Waveforras 

SynfeesizedTe^ 
Point  Wavef<Hms 

TA 

01001 

TB 

01001 

LA 

01001 

LB 

01001 

Bandpass 

IFFT 

Filter 

Bank 

Peak 

Amplitude 

0.6500 

e-7 

|B 

0.0500 

e-7 

BH 

Peak 

Derivative 

0.9768 

c+8 

0.7755 

6+8 

7.3808 

e+8 

4.7091 

6+8 

2.3151 

6+8 

1.5225 

6+8 

Peak 

Impulse 

0.1897 

e-9 

0.5921 

e-9 

m 

Rectified 

Impulse 

0.6700 

e-7 

0.0700 

e-7 

0.0200 

e-7 

0.0600 

e-7 

0.0200 

e-7 

Root  Action 
Integral 

0.0887 

0.0725 

0.5496 

0.4017 

0.1594 

0.1272 
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V.  CONCLUSIONS 


In  this  thesis  the  bandpass  IFFT  filter  bank  and  the  multirate  digital  filter  bank 
techniques  were  developed  to  analyze  and  synthesize  the  component  waveforms  to  produce 
a  synthesized  test  point  waveform  that  retains  the  basic  properties  of  the  component 
waveforms.  The  component  waveforms  are  recorded  by  two  separate  measurement 
systems  (HPW  and  CSW)  under  two  different  aircraft  orientations  (parallel  and 
perpendicular).  The  HPW  waveforms  cover  a  (reliable)  frequency  range  of  0.5-100  MHz 
while  the  CSW  data  are  considered  reliable  over  80-1000  MHz.  This  leaves  a  spectral 
overlap  between  the  HPW  and  CSW  waveforms  in  the  range  of  80-100  MHz.  The  HPW 
waveforms  ,  are  recorded  in  the  time  domain  with  a  sampling  period  of 

9 

At  =  0.4674  X  10  sec.  The  CSW  data  are  recorded  in  the  frequency  domain  in  six 
separate  subbands  with  different  bin  widths. 

In  order  to  produce  a  synthesized  test  point  waveform,  the  following  steps  were 
performed.  The  HPW  and  the  CSW  data  were  transformed  to  be  in  the  same  domain  (time 
or  frequency)  using  the  bandpass  IFFT  techniques.  The  sampling  rates  and  bin  widths 
were  converted  by  inteipolation/decimation  techniques  to  force  them  to  be  the  same  for  all 
waveforms  and  subbands.  Data  from  two  orientations  were  then  combined  using  the  SMA 
algorithm.  Finally,  the  spectrally  overlapped  subbands  were  combined  using  spectral 
concatenation  and  symmetric  overlap  approaches. 

The  tree-structured  filter  bank  method  of  wideband  signal  analysis  provides  an 
efficient  method  by  dividing  a  wideband  signal  into  smaller,  more  manageable  subbands. 
Care  should  be  taken  to  ensure  that  sufficient  number  of  data  points  are  available  for  a  given 
filter  order.  Insufficient  data  length  may  lead  to  end  effects  in  filtering.  In  designing 
filters,  effort  should  be  made  to  minimize  side-lobe  and  “aliasing”  effects.  No  technique 
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is  available  to  completely  eliminate  the  loss  of  information  due  to  aliasing  effect,  but  it  can 
be  minimized. 

Since  the  true  test  point  waveform  could  not  be  physically  measured,  it  was  not 
possible  to  determine  the  accuracy  of  the  synthetic  test  waveforms  produced  in  this  work. 
"N^sual  observations  and  norm  attribute  measures  were  used  to  evaluate  the  closeness  of  the 
synthesized  test  point  waveform  to  its  component  waveforms.  The  results  obtained  from 
the  bandpass  IFFT  and  multirate  filter  bank  methods  were  identical  when  visually 
compared.  The  norm  attribute  values  provide  a  means  for  numerical  comparison,  and 
these  values  for  the  synthesized  test  point  waveform  were  generally  close  to  the  maximum 
norm  attribute  values  of  the  component  waveforms. 
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